Astronomy & Astrophysics manuscript no. shortQUBIC 


©ESO 2010 


October 5, 2010 





QUBIC: The QU Bolometric Interferometer for Cosmology 

The QUBIC collaboration, E. Battistelli 5 , A. Bau 6 , D. Bennett 12 L. Berge 3 , J.-Ph. Bernard 2 , P. de Bernardis 5 , 
A. Bounab 2 , E. Breelle 1 , E.F. Bunn 10 , M. Calvo 5 , R. Charlassier 1 , S. Collin 3 , A. Cruciani 5 , G. Curran 12 , 
L. Dumoulin 3 ,A. Gault 9 , M. Gervasi 6 , A. Ghribi\M. Giard 2 , C. Giordano 5 , Y. Giraud-Heraud 1 , M. Gradziel 12 , 
L. Guglielmi 1 , J.-Ch. Hamilton 1 , V. Haynes 7 , J. Kaplan 1 , A. Korotkov 8 , J. Lande 2 , B. Maffei 7 , M. Maiello 13 , 
S. Malu 11 , S. Marnieros 3 , S. Masi 5 , A. Murphy 12 , F. Nati 5 , C. O'Sullivan 12 , F. Pajot 4 , A. Passerini 6 , S. Peterzen 5 , 
F. Piacentini 5 , M. Piat 1 , L. Piccirillo 7 , G. Pisano 7 , G. Polenta 51415 , D. Prele 1 , D. Romano 5 , C. Rosset 1 , M. Salatino 5 , 
A. Schillaci 5 , G. Sironi 6 , R. Sordini 5 , S. Spinelli 6 , A. Tartari 6 , P. Timbie 9 , G. Tucker 8 , L Vibert 4 , F. Voisin 1 , 

R.A. Watson 7 , and M. Zannoni 6 

1 APC, Universite Paris Diderot-Paris 7, CNRS/IN2P3, CEA, Observatoire de Paris, 10, rue A. Domon & L. Duquet, Paris, France. 

2 Centre d' Etude Spatiale des Rayonnements, CNRS/Universite de Toulouse, 9 Avenue du colonel Roche, BP 44346, 3 1028 Toulouse 
Cedex 04, France. 

3 Centre de Spectroscopic Nucleaire et de Spectroscopic de Masse, UMR8609 IN2P3-CNRS, Univ. Paris Sud, bat 108, 91405 Orsay 
Campus, France. 

4 Institut d'Astrophysique Spatiale, Universite Paris-Sud, Orsay, 91405, France 

5 Dipartimento di Fisica, Universita di Roma "La Sapienza", Roma, Italy. 

6 Dip. di Fisica "G.Occhialini" - Universita degli Studi di Milano-Bicocca Piazza della Scienza, 3 - 20126 Milano, Italy. 

7 School of Physics and Astronomy, The University of Manchester, Alan Turing Building, Oxford Road, Manchester Ml 3 9PL, UK. 

8 Brown University, Providence, RI 02912, USA. 

9 University of Wisconsin-Madison, Madison WI 53706, USA. 

10 Physics Department, University of Richmond; Richmond, VA 23173, USA. 

11 Raman Research Institute, Sadashivanagar, Bangalore 560 080, India. 

12 Department of Experimental Physics, National University of Ireland Maynooth, Maynooth, Co. Kildare, Ireland. 

13 Universita degli Studi di Siena - Rettorato, Via Banchi di Sotto 55, 53100 Siena ITALY 

14 ASI Science Data Center, c/o ESRIN, via G. Galilei, 1-00044, Frascati, Italy 

15 INAF,Osservatorio Astronomico di Roma, via di Frascati 33, 1-00040, Monte Porzio Catone, Italy 

Received ; accepted 

ABSTRACT 

Context. One of the major challenges of modern cosmology is the detection of B-mode polarization anisotropics in the Cosmic 
Microwave Background. These originate from tensor fluctuations of the metric produced during the inflationary phase. Their detection 
would therefore constitute a major step towards understanding the primordial Universe. The expected level of these anisotropies is 
however so small that it requires a new generation of instruments with high sensitivity and extremely good control of systematic 
effects. 

Aims. We propose the QUBIC instrument based on the novel concept of bolometric interferometry, bringing together the sensitivity 
advantages of bolometric detectors with the systematics effects advantages of interferometry. 

Methods. The instrument will directly observe the sky through an array of entry horns whose signals will be combined together 
using an optical combiner. The whole set-up is located inside a cryostat. Polarization modulation will be achieved using a rotating 
half-wave plate and the images of the interference fringes will be formed on two focal planes (separated by a polarizing grid) tiled 
with bolometers. 

Results. We show that QUBIC can be considered as a synthetic imager, exactly similar to a usual imager but with a synthesized beam 
formed by the array of entry horns. Scanning the sky provides an additional modulation of the signal and improve the sky coverage 
shape. The usual techniques of map-making and power spectrum estimation can then be applied. We show that the sensitivity of 
such an instrument is comparable with that of an imager with the same number of horns. We anticipate a low level of beam-related 
systematics thanks to the fact that the synthesized beam is determined by the location of the primary horns. Other systematics should 
be under good control thanks to an autocalibration technique, specific to our concept, that will permit the accurate determination of 
most of the systematics parameters. 

Key words. Cosmology - Cosmic Microwave Background - Inflation - Instrumentation - Bolometric Interferometry 



1. Introduction 

This article describes the proposed QUBIC experiment, a 
Bolometric Interferometer designed to put tight constraints 
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on the Cosmic Microwave Background B-mode polarization 
anisotropies. These odd parity polarization anisotropies are gen- 
erated by primordial gravitational waves (and by lensing of even 
parity polarization at small scales). Detection of these waves 
would represent a major step towards understanding the infla- 
tionary epoch that is believed to have occurred in the early 
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Universe. Tensor modes (primordial gravitational waves) in the 
metric perturbation are indeed a specific prediction of infla- 
tion. The measurement of the corresponding B-mode polariza- 
tion anisotropics would therefore be a smoking gun for inflation. 
A detection would reveal the inflationary energy scale, which 
is directly related to the amplitude of this signal. The tensor 
to scalar ratio r is however expected to be small (smaller than 
0.2 from today's best indirect measurement - the contribution 
of the tensor modes to the temperature and E-mode polariza- 
tion anisotropy) so that the quest for the B -modes is a major 
experimental challenge. Such a small signal (at least an order of 
magnitude below the E-mode signal) justifies the new generation 
of instruments operating from the ground or from balloons (be- 
fore a potential dedicated satellite mission) with unprecedented 
sensitivity and control of systematics. From this perspective, we 
propose the QUBIC experiment, making use of the novel tech- 
nique of Bolometric Interferometry, bringing together the advan- 
tages of bolometric detectors in terms of sensitivity (availability 
of large arrays of background-limited detectors with wide band- 
width) and of interferometry in terms of control of systematic 
effects (clean optics with low induced polarization, low side- 
lobes and therefore low ground pickup, well defined synthesized 
beam, small impact of individual primary beam differences) with 
a large number of detectors which can be replicated quite simply. 

QUBIC will be composed of interferometer modules oper- 
ating at three different frequencies (97, 150 and 220 GHz) with 
25% bandwidth. Each module will respectively comprise 144, 
400 and 625 receiver horns whose signals will be correlated to- 
gether using an optical combiner located inside each module's 
cryostat. After splitting using a polarizing grid, the interference 
fringes will be imaged with two 900 element bolometer arrays 
cooled to 100 mK in order to achieve background-limited sensi- 
tivity. The use of a rotating half- wave plate and polarized focal 
planes allows one to directly reconstruct the "synthetic images" 
of the /, Q and U Stokes parameters observed through the pri- 
mary beam. The TT, TE, EE and BB power spectra can be re- 
constructed from the Stokes parameters' synthetic images using 
standard techniques. The usual calibration and a novel autocali- 
bration technique specific to bolometric interferometry (making 
use of the redundancy of the array providing multiple replica- 
tions of the same baselines) will allow QUBIC to achieve un- 
precedented control of systematics along with a sensitivity com- 
parable to that of more traditional imaging polarimeters. 

Interferometers have a long history of successful measure- 
ments of the Cosmic Microwave Background. In fact, inter- 
ferometers were the first instruments to detect the polarized 
anisotropies flKovac et al., 2002| and [Readhead et al., 2002] ). 
They are generally recognized to offer certain advantages with 
respect to imaging systems for controlling systematic effects (in- 
sensitivity to 1/f noise, clean optics). They also offer a straight- 
forward way to access the angular power spectrum of the sig- 
nal on the sky as they provide a direct measurement of the 
Fourier modes of the sky observed through the primary beam. 
Coherent interferometers require receivers that use either ampli- 
fiers or, for higher frequencies, SIS mixers used in a heterodyne 
arrangement. Any coherent receiver adds an intrinsic and irre- 
ducible amount of noise to the observed signal, preventing het- 
erodyne interferometers from reaching background-limited sen- 
sitivity. Such interferometers are hard to scale to a large number 
of receiving antennas because of the complexity of the correla- 
tors, which must measure the correlations from all possible pairs 
of antennas. The limited bandwidth of these systems and the re- 
quirement to measure polarization increases the complexity. 



All these reasons have led a number of teams to choose 
imaging instruments rather than heterodyne interferometers for 
the next generation of CMB polarimeters. Bolometers are nat- 
urally wide band detectors that have negligible intrinsic noise 
when cooled to sub-Kelvin temperatures. They can be designed 
to operate at the background-limit from the ground or space. 
Up to now they have only been used with imagers which can 
now populate their focal planes with bolometer arrays with up 
to thousands of pixels. Imaging is a fully mature technique, 
widely used to observe Cosmic Microwa ve Background tem- 
perature and polarization anisotropies (Reichardt et al., 2009, 
|The QUAD collaboration,~2009l |Takahashi et al., 20081 ). There 
are however a number of sources of systematic effects that are 
hard to control with imagers (ground pickup and time vary- 
ing systematics such as atmospheric contamination for instance) 
and may prevent such instruments from ultimately achieving the 
exquisite sensitivity required for the detection of the primor- 
dial B -modes. Combining the advantages of interferometry in 
terms of systematics and direct observation of the spatial cor- 
relations of the sky on the one hand with those of imaging in 
terms of sensitivity in the other would obviously be of great in- 
terest for the next generation of Cosmic Microwave Background 
instruments (including future satellite missions). A wide band, 
high frequency and background- limited bolometric interferome- 
ter that could be scaled to a large number of receiving elements 
using large bolometers arrays would potentially be an excellent 
candidate for detecting such a small signal as the primordial B- 
modes. We describe in the present document the design we pro- 
pose for such an instrument. 



2. Bolometric interferometer concept 

The bolometric interferometer we propose with the QUBIC in- 
strument is the millimetric equivalent of the first interferometer 
dedicated to astronomy: the Fizeau interferometer. It was ob- 
tained by placing a mask with two holes at the entrance pupil of 
a telescope. Fringes were then observed at the focal plane of the 
telescope. In our case (see Fig. [T] for a sketch) we use an array 
of back-to-back horns acting as diffractive pupils just behind the 
window of a cryostat. The electric field coming from a given sky 
direction experiences phase differences due to the distance be- 
tween the input horns. The back horns re-emit the electric field 
preserving this phase difference inside the cryostat. The inter- 
ference fringe patterns arising from all pairs of horns are then 
formed on the focal plane of an optical beam combiner which 
is actually just a telescope that superimposes the electric fields 
from all the horns at each point of the focal plane. The polariza- 
tion of the incoming field is modulated using a half- wave plate 
located after the back horns. A polarizing grid separates both po- 
larizations towards two different focal planes each measuring a 
linear combination of /, Q and U Stokes parameters. 

The fact that the electric fields from all horns are added 
and then squared and averaged in time using the bolometers 
makes our instrument an adding interferometer (Fizeau com- 
bination), in contrast to radio-interferometers, where the sig- 
nals (visibilities) are obtained by multiplying the electric fields 
from pairs of receivers (Michelson combination) using an ana- 
log or digital correlator. We do not need a large number of com- 
plex and expensive correlators as the correlation between chan- 
nels is naturally achieved with the bolometers: (\E\ + £2| 2 ) = 

(|£i| 2 ) + (|£ 2 | 2 ) + 2 (Re(£x£*)). The first sum, (|^| 2 ) + (|£ 2 | 2 ), 
which is just the power in the primary beam, is just a background 
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Fig. 1. Sketch of the QUBIC concept 



poweiQ The last term, the interference term, is proportional to 
the visibility. This simple beam combiner allows the instrument 
to be scalable to a large number of input horns. 

We plan to install a first QUBIC module at the Franco-Italian 
Concordia station in Dome C, Antarctica within two years. The 
first module will consist of an array of 400 horns operating at 
150 GHz with 25% bandwidth and 14 degree (FWHM) primary 
beams. The optical combiner will have a focal length of ~ 30 cm 
and each of the two focal planes will be comprised of arrays of 
30x30 bare TES bolometers of 3 mm size. The full instrument 
will include modules at three different frequencies (90, 150, 220 
GHz) and will constrain a tensor to scalar ratio of 0.01 in one 
year of data taking at the 90% confidence level. 



to those required for other CMB polarization experiments based 
on the traditional imager concept. 



3.1. Cryostat and cryogenics 

The first cryogenic stage of the cryostat will be provided by a 
pulse-tube cooler that allows continuous operation of the instru- 
ment at a temperature below 4K. It will cool the filters, horns, 
half-wave plate, optical combiner and polarizing grid. The sec- 
ond cryogenic stage supporting the detectors at lOOmK will be 
provided by a dilution refrigerator. The large (~40 cm) window 
on top of the cryostat will use Zotefoam to provide good me- 
chanical strength and low emissivity. 



3. Instrument sub-systems 

We describe in this section the various sub-systems of the 
QUBIC instrument. Most of them are already available within 
the collaboration but a few still require a reasonable amount 
of R&D (ongoing within the collaboration) to be fully opera- 
tional. Note that all components required for QUBIC are similar 



1 Note that this DC term will vary with time if we scan the sky, but 
on very large timescales due the wide primary beam. 



3.2. Cold optics 
Horns 

The requirements for the front-end (primary) horn antennas 
of QUBIC are similar to those of previous and current CMB 
imaging experiments. The antennas should have low return loss 
and have well-understood beams characterized by low cross- 
polarization, and low sidelobe levels. They should have good far- 
field beam circularity and operate over bandwidths up to about 
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30%. Because of these stringent optical performance require- 
ments, single-mode corrugated horns have to date been favored 
in CMB experiments. This type of horn has already been shown 
to meet the exacting requirements listed above and is readily 
manufacturable at the QUBIC observing frequencies. Single- 
mode horns allow only one orthogonal pair of coherent fields 
to propagate and so have well-understood polarization charac- 
teristics. Few-moded horns, while they have the advantage of in- 
creasing throughput (Murphy et al., 2010] ), are likely to scramble 
polarization information and may not meet the beam circularity 
or sidelobe criteria so easily. 

The length of the horns (if, for example, weight needs to be 
minimized) and the sidelobe levels can be reduced by adjust- 
ing the horn profile away from the stand ard conical shape (e.g. 
IQ'Sullivan et al., 2008l|Maffei et al., 2000| . Shaping can also be 
used to increase the aperture efficiency of a horn by moving the 
phase centre close to the aperture (as could a phase-flattening 
lens). Careful design of waveguide transitions keeps the return 
loss very low. We are considering either of the two types of 
horns developed for the Clover project. At 150 GHz "Ultra- 
Gaussian" horns (so- called because 99.9% of power is in the 
free-space fundamental Gaussian mode) were designed (Ade et 
al., 2009). In these horns the HE12 mode is deliberately excited 
in a cosine-squared profile section of the horn and then brought 
into phase with the dominant HEn mode in an extended parallel 
front section. At 97 GHz corrugated Winston-like profiled horns 
(Maff ei etal.72 004) were used and shown to produce very low 
spillover (Grimes et al., 2009| ). 

Filters 

Spectral filters will be located inside the dewar to perform sev- 
eral roles. First, the selection of the spectral band of observation 
is achieved by the combination of the back-to-back waveguide 
section (frequency cut-on) and metal-mesh interference filters 
(frequency cut-off). Second, metal mesh interference filters lo- 
cated at the entrance aperture of each thermal stage reduce the 
radiation load on the cryogenic system and the detectors. The ra- 
diation load reaching the detectors is not only a source of stray- 
light but will also lead to a poorer sensitivity if not well under 
control. 

- 300K stage: The dewar window will already act as a filter by 
reflecting and absorbing most of the short wavelength radia- 
tion, up to the near-infrared. 

- Intermediate thermal stage filters in front of the back-to-back 
horns: Each screen aperture will have at least one low-pass 
edge filter with possibly an associated thermal (IR) filter 
to reduce the thermal loading. These filters and the 300K 
stage window and filter, will all have to be carefully designed 
and characterized. Being located in front of the back-to-back 
horn apertures, they can affect the beam shape of the instru- 
ment. Previous measurements have shown that the best com- 
ponents have a minimal impact on the beam shape and an in- 
duced instrumental polarization lower than -35dB. However 
they will have to be accurately characterized, typically down 
to -50dB. In order to clear the main beam of all the horns, 
and taking into account the number of pixels per instrument, 
a clear aperture diameter of 400mm is required for all these 
filters and windows. 

- Filters located after the back-to-back horns: Probably two 
thermal stage apertures will be between the back-to-back 
horns and the detectors (4K and lOOmK). Although the fil- 
ters located at these apertures will be manufactured with the 



same attention as the previous ones, because the polarization 
will have been already modulated, they will have a lesser im- 
pact on systematic effects. The focus will be on achieving a 
high out-of-band rejection in order reduce stray-light. The 
spectral band defining filter will be located on the aperture 
of the box shielding the detectors in order to be as close as 
possible to these. 

Switches 

The switches are required during the calibration procedure but 
are not used during data acquisition. They are used as shutters 
that are operated independently for all channels. The idea is to 
modulate a given visibility measured by a given pair of horns at 
one time in order to compare it with other pairs of horns measur- 
ing the same (equivalent or redundant) baseline. This has proven 
to be an extremely powerful self-calibration tool in interferom- 
etry (see section [6] for details). The switches will metal pins 
or films that close the waveguide section between the back-to- 
back horns. They will be operated using miniature electromag- 
nets commanded from the outside of the cryostat. The switching 
does not need to be fast but will need to be reliable. 

Half-wave plate 

The rotating half-wave plate modulates the signal coming from 
the sky. It is therefore a critical component that needs to be 
designed carefully and tested extensively. The systematic ef- 
fects that could be generated by the h alf-wave plate are re - 
duced when cooled to low temperature ( Salati no et al., 2 010a). 
The inhomogenities of the half-wave plate are however an is- 
sue for such large modulator. Possible solutions to manufac- 
turing a broadband half-wave plate include either birefringent 
plates like sapphire fo r example flPisano et al., 2006| ) or a stack 
of metal mesh filters ( |Pisano et al., 2008| ). The rotating mecha- 
nism will be based on the design made for the PILOT experiment 
( Salati ncTet al., 2010a| ), allowing for either a continuous rotation 
speed or a step rotation with a precision better than 0.1 deg. 

Optical combiner 

Signal combination in QUBIC will be performed by means of 
an optical system that transfers the fields radiated by the pri- 
mary back-to-back horn array to a detector plane where Fizeau 
interference fringes can be observed. The QUBIC system must 
satisfy the following basic requirements: 

- The combiner is an optical system in which rays launched 
at a given angle from the re-emitting (back) horn array are 
focused to a single point on the focal plane. In this way (for 
an ideal optical system) equivalent baselines will produce 
identical fringe patterns. 

- The limit on the total number of bolometers that can be pro- 
duced, together with the sampling requirement for at least 
two bolometers per fringe, constrains the equivalent focal 
length of the combiner to range from 200 to 300 mm. Since 
the back-to-back horn array size is 240 mm in diameter for 
the 150 GHz instrument, the combiner will be a very fast 
system (small F#). 

- We can use mirrors or lenses as optical components. Our 
choice will be between an on-axis lens system or an off-axis 
mirror system that can avoid the shadowing of any baselines 
(on-axis mirrors would result in very high levels of trunca- 
tion). On one hand the behavior of mirrors can be extremely 
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well characterized using physical optics, but on the other an 
off-axis system will introduce aberrations. 
- The system has to be compact enough to be placed in a cryo- 
stat of about lm 3 . 



The collaboration is investigating different layouts for the com- 
biner, possibly including lenses, starting with classical astro- 
nomical telescope designs such as Gregorian and Cassegrain, 
and bearing in mind the constraints above. We find that CATR 
designs (Compact Test Range, crossed Dragone configuration - 
|Tran et al., 20 08), which have proven suitable elsewhere, are not 
quite feasible given our very short focal length and wide field- 
of-view. We are therefore currently studying the optimization of 
an off-axis Gregorian telescope. 



3.3. Detector chain 

The focal plane will be covered with a polarized filled array in 
order to sample properly the image. The simplest way to separate 
polarization is to use a polarizing grid which directs the beam to 
two detector arrays, one in transmission and one in reflection. 
Taking into account the constraints from the optical combiner 
and the size of the secondary horns leads to the need for about 
900 detectors of 3mm size each for the 150 GHz module. The 
main performance requirements for the detectors are their NEP 
and time constant. With a background power of few pW, the de- 
tector NEP needs to be lower than about 5 - 10aW.Hz" 1/2 at 
150GHz. The modulation speed of the half wave plate of a few 
Hz constrains the detector time constant to be shorter than 10ms. 
These requirements are easily reached with lOOmK TESs which 
are the current baseline. 

The detector assembly sub-system will be based on NbSi 
transition edge sens ors (TES) read out with a time domain 
multiplexing scheme ( [Pajot et al., 2008, Pi atet al., 20081 ). While 
other technologies like KIDs could be used, TES array technol- 
ogy is currently more fully developed. The use of NbSi alloy 
as the thermal sensor offers several advantages with respect to 
classical bi-layers: the sensors are more homogeneous on a large 
array, the thermal response is high and the noise properties of 
NbSi alloy is intrinsically better than bilayers. A first design of 
a detector array for the focal plane is shown in Fig. [2] 

The TES readout electronics use SQUIDs as a first amplifier 
stage followed by a 4K SiGe ASIC that contr ols the multiplexing 
and amplifies the signal from the SQUIDs ( Vois in et al., 2 008, 
Prel e et al., 20 09). The advantages of such a system are a simpli- 
fication of the architecture, miniaturization of the readout elec- 
tronics and also immunity to high energy particles along with an 
overall low power consumption and very low noise properties. 
We have demonstrated this technique with a 24 pixel readout in 
a 24:1 multiplexing scheme (done in two steps: 8:1 by SQUIDs 
stage and 3:1 with the ASIC). The current design is based on 
a new ASIC under development that will be able to readout 128 
detectors with a 128: 1 multiplexing factor (32: 1 SQUIDs and 4: 1 
ASIC). The full detector array will be therefore read out with 7 
of these ASICs. 



Mount 

The instrument will be installed on a classical alt-azimuthal 
mount allowing for three axis motion in order to both scan in 
azimuth and elevation and to rotate the instrument around its op- 
tical axis. 




Fig. 2. Mask design for a 128 elements TES array on a 2 inch 
silicon wafer. We have also displayed a zoom of one of the de- 
tectors on the bottom-right. 

4. Data analysis for bolometric interferometry 

4.1. Synthetic image 

One can express the bolometer measurements as a linear combi- 
nation of the visibilities (basically Fourier modes) of the sky ob- 
served through the primary beam of the horns. The coefficients 
of the linear combination of visibilities arriving at each bolome- 
ter are just related to the phase-differences between different 
channels imposed by the optical combiner. This was explored in 
detail in (Charlassie r et al., 20 09) where we have shown that the 
visibility reconstruction can be done in an optimal way. Working 
with the visibilities, although it may seem logical as we are deal- 
ing with an interferometer, may not be the best approach how- 
ever. It is better to use our instrument as a synthetic imager. As 
will be shown in this section, the image in our focal plane is 
exactly the quantity usually known in interferometry as the syn- 
thetic image and is usually obtained through an inverse Fourier 
transform of the visibilities. In our case, the visibilities are not 
the natural observable (in contrast with a multiplicative interfer- 
ometer) as we have immediate access to the synthetic image. We 
can however recover the visibilities through a Fourier transform 
of the synthetic image. The synthesized beam we convolve with 
the sky is determined by the configuration of the primary horns, 
its resolution being basically given by the largest separation be- 
tween horns. Once this synthesized beam is known, one can use 
our instrument as a simple synthetic imager by scanning the sky 
as with an usual imager, making maps and calculating power 
spectra using standard techniques. It is however completely an 
interferometer in the sense that the synthetic beam is well de- 
fined by the location of the primary horns. One of the most im- 
portant advantages of having an interferometer is the fact that, 
using the redundancy of the baselines, one can internally cali- 
brate the instrument accurately, including systematic effects such 
as gains and cross-polarization. This topic will be discussed in 
section |6] 

The signal on the bolometers as a function of time i^j 



R(d p , t) = S iidp) ± cos(AcL>t)S Q (d p ) ± sin(4ajt)S v (d p ) 



(1) 



2 This is however an approximate expression as multiple reflexions 
will modify it, producing terms at oj and 3oj that have to be dealt with 
flSalatino et al., 2010a|rSalatino & de Bernardis, 2010b| ). 
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where the ± is + for one of the focal planes (polarized in one di- 
rection) and - for the other one polarized in the other direction. 
d p is the location of the detector in one of the focal planes and oj 
is the rotation frequency of the half-wave plate. The combination 
of hase- sensitive detection at the half- wave plate modulation fre- 
quency and scanning across the sky allows recovery of the three 
focal plane images, Si(d p ), S Q(d p ) and Su(d p ), corresponding 
to the three Stokes parameters. 

For the X Stokes parameter the synthetic image can be di- 
rectly expressed as a function of the observed sky: 



Sx(d p ) = J X(n)B p s (n)dn 



(2) 



where B p s is the synthesized beam of the interferometer for de- 
tector p that is formed by the arrangement of the array of input 
horns. 



4.2. Synthesized beam 

The synthesized beam is just the sum of all wide band fringe 
patterns formed by all pairs of horns modulated by the primary 
beam on the sky and by the secondary beam on the focal and 
integrated over the bolometer surface: 



B p s (n) = B prim (n) J J B SQC (d) 



^ exp ^27^ • — n 



(3) 



J(v)@(d - d p )dvdd 



where Df is the equivalent focal length of the optical combiner, 
#prim is the primary beam on the sky of the input horns, B sec (d) 
is the beam of the reemitting horns as seen on the focal plane 
(d labels the position on the focal plane). The integration is per- 
formed over the bandwidth of the instrument J(y) (v is the fre- 
quency and X - c/v the wavelength) and over the surface of 
the bolometer indexed by p and modeled with the top-hat like 
function ®(d - d p ). The sum in the integral is performed over 
the input horns with position given by x t . We have assumed in 
the above expression that all horns have the same primary and 
secondary beams so that they can be factorized. A more real- 
istic case can be easily accounted for in the calculation of the 
primary beam, provided the fact that individual beams are well 
known (from calibration on point sources for instance). 

If one had an infinite number of primary horns, the synthetic 
beam would be a Dirac peak in the center of the field of view of 
the primary beam. We have chosen to place the primary horns on 
a compact square array as shown in Fig. [3] in order to achieve a 
maximal redundancy of the baselines (see Fig. [4} defined by each 
possible pair of horns. Such a redundancy was shown to be cru- 
cial to ensure an optimal sensitivity of the bolometric interferom- 
eter (Charlass ier et al., 2009| ) allowing for what we have called a 
"coherent summation of equivalent baselines". This means that 
all equivalent (or redundant) baselines need to experience the 
same phase difference configuration on the focal plane in order 
to maximize the signal to noise ratio of the visibilities, or of the 
synthetic images. Such a condition can be easily understood: on 
the focla plane, each baseline contributes a sinusoidal fringe pat- 
tern with a spatial frequency determined by the baseline length. 
A second baseline equivalent to the first one would contribute 
with the same pattern shifted if both baselines were not "phase- 
locked". With a large number of equivalent baselines N Qq , the 
fringe signal would be washed out as 1 / ^jN Qq while if all base- 
lines are coherent in phase, the fringe signal is just summed N Qq 
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Fig. 3. 20x20 horn array configuration for the 150 GHz QUBIC 
module. Each horn has a beam of 14 degree (FWHM), has an 
internal radius of 5.8 mm, and has 1 mm thick walls. 



150 



100 



50 




-50 



-100 



-150 



::::: ::::: 
••••••••••••• ••••• •••••••••••••• 

•••••••••••• ••••••••• ••• 



76p. different :ba$elin$s. 
79809 pairs- in total- - ■ 



Neq-=-3€9m 



-150 



-100 



-50 



50 



100 



150 



Fig. 4. uv -plane coverage of the array shown in Fig. [3] Each point 
corresponds to a different baseline and is highly redundant in our 
array. The number of equivalent baselines is shown in colors. 
The multipoles can be calculated from the baseline length by 
multiplying by 2n (this is exact in the flat sky approximation, but 
we are not in this range, so the correspondence is approximate). 



times without being washed out. The use of an optical combiner 
automatically ensures that equivalent baselines are summed co- 
herently (Charlas sier et al., 2009| ). A compact square array of 
primary horns ensures a large number of equivalent baselines 
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(hence a good sensitivity) while allowing for an almost maximal 
density of horns within the same aperture window. 

With such a redundant array, in the monochromatic and 
point-like detector limit, the integrals in Eq. [3] disappear and the 
synthetic beam can be analytically calculated: 

B p s (n) = B vrim (n)B SQC (d p ) 



(4) 



where P = 20 is the number of horns on a side of the square 
array, Ax is the distance between two horns and the indices x 
and y denote projection along the respective axis in the focal 
plane reference frame. £ m [ n ~ 2n^f ~ 43 is the minimum mul- 
tipole accessible with this array and a single field observation 
while ^ max — 27r^p ~ 867 is the largest multipole accessi- 
ble (although with a small number of equivalent baselines, and 
hence a very poor sensitivity - the actual maximum for well- 
measured multipoles is around 200). The maximum multipole 
accessible defines the angular resolution of the synthetic beam 
to be ~ — 0.4 degrees while the fact that the minimum mul- 
tipole accessible is £ m [ n implies that this small synthetic beam 
replicates over the sky with peaks separated by ~ j^- ~ 8.5 de- 

£min 

grees. This replication is however significantly reduced by the 
apodization of the primary beam that defines the field of view of 
our interferometer (we have chosen 14 deg FWHM). The peaks 
of the synthetic beam are shifted relative to the center of the 
field of view for different detector locations in the focal plane. 
The primary beam however remains fixed so that the relative in- 
tensity of the synthetic beam peaks changes from one detector to 
another. 

The synthetic beam is shown in Fig.[5]for monochromatic or 
broadband cases and point-like or finite- size detectors (located 
at the center of the focal plane) showing the effect of the integra- 
tion over bandwidth and detector surface. We have used an array 
of 400 primary horns with 14 degree (FWHM) primary beam 
and a combiner with an equivalent focal length of 30 cm, a cen- 
tral frequency of 150 GHz and 25% bandwidth. We also show in 
Fig. [6] horizontal cuts of the maps in Fig. [5] showing the profiles 
of our synthesized beam. We see that integration over the pixel 
size slightly degrades the resolution of the central peak from 0.33 
degree to 0.54 while integration over bandwidth only affects sig- 
nal far from the center of the field of view. We also note that as 
expected, the replicated peaks are separated by ~ 8.5 degrees. 

4.3. Power spectrum from the synthetic image 

In the previous section, we have shown that the direct observable 
of a bolometric interferometer such as QUBIC is the synthetic 
image, usually obtained in traditional interferometry from the in- 
verse Fourier Transform of the visibilities directly measured by 
the set of horns and correlators. Our bolometric interferometer is 
therefore very different from other interferometers in this respect 
and somehow closer to an imager whose instantaneous observ- 
able is also the sky convolved with the instrumental bearrj^] 

The visibilities in an interferometer are quantities that are 
defined as the Fourier transform of the observed field pro- 
jected on the plane tangent to the observing direction. These 
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3 In an imager, this beam is formed by the antenna that couples to 
each detector projected on the sky through the telescope (reducing the 
horn beam from a few degrees to a few arc-minutes depending on the 
telescope diameter). 
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Fig. 5. Gnomonic projections of Healpix ( [Gorski et al, 2005| ) 
maps of the QUBIC synthesized beam for a monochromatic in- 
strument with point-like (top-left) or 3 mm detectors (top-right), 
a 25% bandwidth instrument with point-like (bottom-left) and 3 
mm detectors (bottom-right). All images were obtained with a 
compact square 400 horns array with 14 degrees FWHM beams 
at 150 GHz. The effect of the integration over bandwidth and de- 
tectors can easily be seen, the former only affects the shape far 
from the center while the latter slightly enlarges the peaks. The 
black circles are at 5 and 10 degrees from the center. 



are just the modes on the sky if the primary is small enough 
so that the flat sky approximation is valid. With a 14 degree 
primary field of view, this approximation is clearly not valid. 
It is therefore not straightforward to combine visibilities taken 
while observing different directions of the sky with a wide field 
of view interferometer ( although it has been explored in detail 
in |Bunn et White, 20 06). With our synthetic imager approach, 
the problem is naturally solved: as the synthetic image is simply 
the convolution of the sky through the synthetic beam, one can 
directly relate the covariance matrix of the bolometer data to the 
angular power spectrum. Let's introduce /3{ m (p, no, 6q): 



B p s (n - no, 6 )Y £m (n)dn 



(5) 



where B p (n-no, 0q) is the synthesized beam for detector number 
p in the focal plane when the intrument points towards no with a 
pitch angle 6q (rotation around the no direction Q 



4 Note that /3 £m (p, #o) is not the spherical harmonic coefficient of 
B p s (n - no, #o) (which would have required a Y* m ), but the conjugate of 
the spherical harmonics transform of B p s *(w - w , 6q). 
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Fig. 6. Horizonthal cut of the synthesized beam maps shown in 
Fig. [5] The FWHM of the central peak is obtained through a 
Gaussian fit to the profile. One sees the slight loss of resolu- 
tion due to the finite detector size while the effect of bandwidth 
(dashed lines) can be seen on the replicated peaks, far from the 
center of the field of view. The next peak not displayed in the 
figure (at ~ 17 degrees from the center) peaks at ~ 0.015 for 
the blue and green lines and ~ 0.006 for the orange and red 
lines. The following one (at an angle of -25.5 deg) peaks around 
~ 5 x 10" 5 . These numbers assume a perfect gaussian primary 
beam and might be larger in the presence of sidelobes in the pri- 
mary beams. 

It is straightforward to show that the total intensity synthetic 
image (measured with bolometer p when pointing towards no 
with pitch angle 0q) is: 



S /(/?, no, G ) = ^ a tn fi tm {p, w , O ) 



(6) 



Cm 



where the a^ m are the underlying sky spherical harmonic expan- 
sion coefficients. The covariance matrix of the bolometer signals 
can be explicitly written as the following: 

(s fanned 'S*(p\n u 0j)) = Y J CiW i (p,p\n k ,n h 6 u e j ) (7) 

i 

where we have introduced the window function: 

W € (p, p\ n k , n u U Oj) = J^fo **> 9 i)fiL(P'> n i> e i> ( 8 ) 

m 

Similar window functions can be calculated in a straightforward 
manner for polarization of course, so that the covariance matrix 
of the synthetic images of /, Q and U can be expressed as a func- 
tion of these window functions and the power spectra C rr , C TE , 
C EE , C BB (as well as C TB and C EB generally required for consis- 
tency checks). If one observes a limited number of fields with a 
limited number of pitch angles so that the number of bolometer 
data samples involved is not too large, the power spectrum can 
be estimated using a brute force maximum likelihood using the 
window functions in Eq.[8] 

It is useful to have a look at the window functions as they 
give a good feeling of the sensitivity of the instrument to various 
multipoles. The diagonal parts of the window functions (p = p\ 
Hk = nu 0i = 0j) are shown in Fig. U\ in the case of either point- 
like or 3 mm detectors, with or without bandwidth smearing in- 
cluded. The peaks that are clearly seen in this figure show that 



Fig. 7. Diagonal (p = p\ % = nu 0i - 0j) window functions 
for QUBIC (in arbitrary units) in the cases of point-like or 3 mm 
size bolometers, neglecting and accounting for 25% bandwidth 
smearing (dashed lines). The two bottom curves (dotted lines) 
show the evolution of the signal-to-noise ratio expected on the 
power spectrum accounting also for the noise penalty coming 
from bandwidth smearing (see Eq. 10). Note that above £ ~ 200, 
the statistical noise becomes dominant although it does not ap- 
pear on this plot. 



we have access to some discrete sample of all the multipoles with 
a resolution in multipole space given by the inverse of the pri- 
mary beam ~ 1 /cr pr i mary . One can actually exactly recover the 
monochromatic, point-like detector window function just from 
the number of equivalent baselines at each multipole and this 
multipole space resolution. One clearly observes the suppression 
of power at high multipoles due to integrating the fringes smaller 
than the size of the detectors as well as the reduction of power 
and smoothing in multipole space due to bandwidth smearing. 
The latter is in perfect agreement with analytical calculations 
performed in visibility space in ([Charlassier et al., 2010| ). One 
should notice that the bandwidth smearing also causes an addi- 
tional signal-to-noise reduction (enhancement of the noise due 
to the reduction of the effective primary beam, see Eq. [TO} 



that is not visible in the window function which just shows the 
amount of signal that is filtered by the instrument. This is de- 
tailed in (Charlassie r et al., 20 10). For this reason we show in 
Fig.[7]the two bottom curves (dotted lines) which display the ac- 
tual evolution of the signal-to-noise ratio as a function of i due 
to the full bandwidth effect. 



4.4. Scanning the sky 

Using QUBIC just as an imager is therefore completely equiv- 
alent to using it as an interferometer. It actually brings in many 
advantages through the possibility of scanning the sky, like usual 
imagers and unlike usual interferometers. Scanning the sky ac- 
tually offers an immense advantage with respect to observing 
single fields thanks to the extra-modulation of the sky signal it 
allows. The half- wave plate polarization modulation shown in 
Eq. [T] does not allow recovery of the total power signal, unless 
one scans the sky, as total power is not modulated by the half- 
wave plate. Scanning the sky is therefore the only way of re- 
covering / in addition to Q and U. Furthermore (and even more 
important), spurious signals coming from reflections of unpo- 
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larized background light signals on the various components in 
the optical path, polarized by the polarizing grid, modulated and 
re-reflected by the half- wave plate, will induce significant distur- 
bances precisely at four times the modulation frequency of the 
half-wave plate, just where the Q and U sky signals are mod- 
ulated flSalatino et al., 2010a| |Salatino & de Bernardis, 2010b| ). 
This m ajor issue with hal f- wave plates has been addressed in de- 
tail in (Joh nson et al., 2007] ) and the removal of this waveplate- 
synchronous disturbance proved to be only possible through 
scanning the sky: this scanning slightly shifted the signal modu- 
lation frequency sideband from the central 4a>, where the distur- 
bance was dominant. In our case, scanning is therefore necessary 
to extract the synthetic images modulated by the half- wave plate. 

4.5. Scanning and Map- making 

Exactly in the same way as for a usual imager, we can scan 
the sky with our synthesized beam with any scanning strategy. 
Each individual measurement results in the convolution of the 
sky through our synthesized beam. The Time-Ordered Data of 
each individual bolometer can then be reprojected on the sky us- 
ing standard map-making techniques. The usual pseudo-power 
spectrum estimation can then be applied to the resulting maps. 
Such simulations are currently ongoing within the QUBIC col- 
laboration. 

We plan to use a scanning strategy based on azimuth scans 
with AAz ~ 10 degrees slowly varying the elevation (ranging 
from 45 to 65 degrees) after a number of scans in order to have 
the center of the field of view scanning an approximate circle 
on the sky. The scanning parameters must be adapted to the 
instrument design. The primary beam of QUBIC is around 14 
degrees; the azimuth scans should therefore be less extended 
than for traditional imagers, to avoid galactic contamination. 
We have chosen a typical azimuth scan of 10 degrees, close 
to the distance between main and secondary peaks of the syn- 
thetic beam. The scanning speed must be fast enough to modu- 
late the / sky signal above the 1 //-noise: assuming the detectors 
are stable over a period of 10 s, it requires the scanning speed 
to be of the order of 2 degrees per second. Finally, the half- 
wave plate should rotate fast enough to allow measurement of 
Q and U before the main peak of the synthesized beam moves 
more than a fraction of a beamwidth. Given the scanning speed 
and the main peak size (around a degree), a rotation frequency 
of a few Hertz is enough. These parameters are actually not 
very far from the ones chosen for previous imaging experiments 
(Takahashi et al., 2008) | Johnson et al., 2007] ), and so should be 
technically achievable. As a result the final field of view (shown 
in Fig. [8]) defined by the sum of the primary beams along this 
circle will be significantly flattened and would achieve a value 
of T] = 1.6 (see Fig.[9]for an explanation of n). 



5. Sensitivity comparison with an imager 

It is of course critical to estimate the sensitivity of an instru- 
ment like QUBIC and to compare it to that of an imager with 
the same number of horns and similar angular resolution. This 
study has b een performed in deta il, neglecting the effect of band- 
width, in (Hamil ton et al, 2 008). A full account of bandwidth 
sm earing and its effect on CMB measurements can be found 
in ( |Charlassier et al., 2010| ) so that one can obtain a general for- 
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Fig. 8. Field of view, or sky coverage, obtained with QUBIC by 
scanning ~ 10 degrees in azimuth with slowly varying elevation 
(from 45 to 65 degrees). The fraction of sky covered is ~ 2% of 
the whole sky with a rather flat profile achieving r] = 1.6. 

mula for the sensitivity of a bolometric interferometer on the 
CMB polarization (E or B) power spectrum: 

BI / 2*i(fl ( 2^NET 2 Q _ x A 

where t is the integration time, £1 = 2^cr^ eam is the solid angle 
subtended by the primary beam (therefore the field of view of 
the instrument), / S k y = £l/2n is the fraction of the sky covered 
by the instrument. Nh is the number of horns installed in a com- 
pact square array, N eq (£) is the number of equivalent baselines 
at multipole NET is the noise equivalent temperature of the 
detectors that scales oc 1 / VAv. ?] in the above equation stands 
for the apodization factor of the field of view n = . For 

a Gaussian primary beam it is equal to 2 but would be one if the 
field of view had a rectangular shape. By scanning the sky with 
the bolometric interferometer, we will achieve a value for r] that 
will be between 1 (optimistic) and 2 (pessimistic). ebi is the over- 
all optical efficiency of the bolometric interferometer. w p t x {i) 
represents the loss in signal to noise due to integration over the 
detector area (this is the ratio of the blue and green curves in 
Fig. [7]). K\(£) is a quantity that defines the effect of bandwidth 
smearing in interferometry (see Charlassie r et al., 2010| for de- 
tails), for a Gaussian primary beam and a Gaussian bandwidth 
with Ay = <x y V27T, it can be analytically calculated: 

*!(/) = S + ^ffi (10) 

where 07 is the resolution of the interferometer in multipole 
space, o~t = 1/cTprimary. Equation [9] is the equivalent for bolo- 
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metric interferometry of the well known Knox formula for im- 
agers flKnox, 1 997): 



AC* m = 



! (2^+l)/ sky A, 



27?NET 2 ^ ! 



Im 



(ID 



where we assume a uniform sky coverage and the beam window 
function is given by = exp + l^beam]- ^ ne P arame ~ 

ter 7] corresponds to the sky coverage uniformity (rj = 1 for a 
uniform sky coverage). 

One can therefore calculate the ratio of the noise parts of 
these equations to have a feeling for the relative sensitivities be- 
tween an imager and a bolometric interferometer. We will as- 
sume that they have the same number of horns and cover the 
same sky fraction (in one single field for the bolometric interfer- 
ometer, and by co-adding maps for the imager). The ratio is: 



AC* m 



ACf 



TJln 



xw piK (£)XKf 2 (e)x 



NET 



NET*, 



Im x^(12) 



We show the result of the above ratio in Fig. [^assuming the same 
NET and optical efficiencies for both instruments. We show 
different cases for the sky apodizations corresponding to ideal 
or more realistic scanning strategies for both the imagers and 
QUBIC. The curve that is the most likely to be realistic is the 
solid-red one corresponding to a realistic sky apodization (qi m = 
1.4 for the imager, taken from QUAD - |Brown, priv. com m. ^| 
and to the scanning strategy presented in Fig. [8| for QUBIC with 
t/bi = 1.6. We observe that the bolometric interferometer seems 
slightly less sensitive than the imager, but the difference is small 
considering the uncertainties in all the factors used for the com- 
parison. So one can conclude that there is no clear advantage 
from the statistical sensitivity point of view to any of the two. 



6. Systematic effects 

The expected weakness of the B-mode signal compared to tem- 
perature and E-mode anisotropies implies that instruments ded- 
icated to this search not only need to be highly sensitive from a 
statistical point of view, but they also need to have an exquisite 
handle on systematic effects that usually cause total power and 
polarization leakage. This reason is the driver for our efforts 
to develop a bolometric interferometer as interferometers are 
known to allow for a very accurate accounting of systematic ef- 
fects. 



6.1. Different systematics in interferometry and imaging 

A first point regarding systematic effects is that we anticipate 
them to be significantly different in our interferometer as com- 
pared to all other experiments, most of which are imagers with 
bolometric detectors. It has been recognized for a long time that 
systematics effects behave in a significantly different manner 
with an interferometer than with an imager ( |Bunn, 2006| ). A non- 
exhaustive list of differences relevant to our case follows: 

- Pointing mismatch and gain differences between detectors 
that are differenced in order to measure Q and U will not be 



5 The value of r/i m = 1.4 quoted for QUAD corresponds to the result 
at the end of the data analysis process and therefore includes all filtering. 
It is hard to compare fairly with the number we quote for QUBIC that 
is an anticipation. 
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Fig. 9. Comparison of the sensitivity to the CMB power spec- 
trum of an imager with 1 degree (FWHM) resolution with a 
bolometric interferometer with a primary beam of 14 degrees 
(FWHM). Each system has 400 horns, 25% bandwidth, and a 
center frequency of 150 GHz. The interferometer has detectors 
of size 3 mm. The dashed curves correspond to an imager with 
a perfectly uniform sky coverage (r]i m = 1) while the solid ones 
are for a more realistic sky apodization including actual scan- 
ning strategy and filtering with r]i m = 1.4 (taken from QUAD - 
|Brown, priv. commT] ). The colors corresponds to various assump- 
tions regarding the bolometric interferometer scanning strategy: 
single field (no scanning) with /7 B i = 2 is in blue and an unre- 
alistic uniform sky coverage (t/ B i = 1) is in green. The present 
scanning strategy for QUBIC corresponds to tjbi = 1.6 is shown 
in red. 



an issue in our case thanks to the use of the half-wave plate 
that allows a measurement of all Stokes parameters with one 
total power detectoi^] 

Cross-polarization coming from the telescope in an imager 
occurs before the half- wave plate and is therefore modulated 
as the sky signal. In our case, the half-wave plate is located 
right after the horns and there is no telescope before. We ex- 
pect the cross-polarization from the window, filters and horns 
to be very small. The cross-polarization induced by the tele- 
scope or the polarizing grid will not be modulated as the sky 
signal and will therefore not affect the signal in a dangerous 
manner. 

The interferometer's scanning strategy can be slow as it does 
not need to scan a large area of the sky (the primary beam 
is already large) so long bolometer time constants (that are 
often far from pure exponentials) will have a lesser impact 
than for an imager. 

Ground pickup from sidelobes induced by the telescope have 
been a strong limitation for imagers from the ground. The 
interferometer primary beams are determined by horns that 
have small intrinsic sidelobes. 

The interferometer's resolution (the synthesized beam) is 
completely defined by the primary horn array and differ- 
ences in the primary beams will have a very small impact 
on the reconstructed signal. The synthesized beam can be 
calculated in several ways (see below) which can be com- 
pared and checked for consistency. The strong limitation in 



6 Note that most upcoming imagers dedicated to B-mode search will 
also use a half- wave plate to modulate the polarization. 
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high-multipoles measurements for imagers due to beam un- 
certainties will therefore be mitigated in our case thanks to 
the extreme accuracy of our knowledge of the synthesized 
beam. 

Having an instrument affected by systematics in a different way 
as the others seems to be a very important step towards having a 
convincing detection of the B-mode signal. 

6.2. Autocalibration 

The main reason to think that systematics might be more con- 
trolable in a bolometric interferometer than in a traditional im- 
ager is, however, different. It is related to a specific autocal- 
ibratio n technique that has be en developed in our c ollabora- 
tion ([Charlassie r et al., in prep. | and Charlas sTer , 2010| ). The au- 
tocalibration allows the determination of most of the system- 
atic effects (as modeled using the Jones matrix formalism for 
instance) independently for all channels from observations of a 
polarized source (whose polarization does not need to be known) 
using the switches that can be seen in the schematic view of the 
instrument (Fig.[T]) between the back-to-back horns. 

This autocalibration technique is inspired by traditional 
interferometry (especially long baseline and optical) where 
the phase is often lost due to atmospheric turbulence 
(see |Pearson and Readhead, 1984] for a review). Two main 
techniques have been used to solve this problem in the 
past. The first technique is based on clo sure phase (see 
for instance Readhead and Wilkinso n, 1978] ) where the un- 
known phases are iteratively reconstructed by forming quan- 
tities where the unknown phases are nulled (for instance 
the product of the three visibilities one can form with 
three receivers). The second, more powerful technique is 
based on the redundancy of the receiver array ([Yang, 1988, 



0.5 



Wieringa, 1991 Noordam and de Bruyn, 1982]). It is rarely 



used because most interferometers have low redundancy as 
they seek high angular resolution and good image recon- 
struction (dense uv-plane coverage) rather than sensitivity. 
The "omniscope" proposed in (Tegmark and Zaldarriaga, 2009a 
and [Tegmark and Zaldarriaga, 2009b ) however relies on redun- 
dancy in the same way we do (it is actually a numerical 
equivalent of our bolometric interferometer, but focused on 
21c m observations) and uses them as a source of autocalibra- 
tion ( [Liu et al, 2010| ) although they are not sensitive to polariza- 
tion. 

The autocalibration technique uses the fact that in the ab- 
sence of systematic effects, equivalent baselines of the inter- 
ferometer should measure exactly the same quantity. Using the 
switches located between the back-to-back horns one can mod- 
ulate on/off a single pair of horns while leaving all the others 
open in order to access the visibility measured by this pair of 
horns alone. By repeating this with a subset of all available base- 
lines (equivalent and different), one can construct a system of 
equations whose unknowns are the systematic effects parameters 
for each channel (as modeled using Jones matrices for instance) 
meaning two polarizations for each primary horn on each of the 
bolometers of the focal plane. One can show that for a large 
enough array of primary horns (at least ~ 20 horns) the system is 
over-constrained and can be solved. No information is required 
on the actual polarization of the observed source except that it 
needs to be polarized. The autocalibration requires the knowl- 
edge of the individual primary beams of each horn that can be 
obtained through scanning an external unpolarized source. Once 
all the couplings and gains are known for each channel and for 
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Fig. 10. Simulation of the autocalibration technique: we show in 
black the set of visibilities used as an input (the x-axis is just 
the label of the different visibilities) generated with a 144 horns 
bolometric interferometer. The dashed red line shows the visibil- 
ities reconstructed without accounting for systematics. The cor- 
ruption is obvious and is completely solved after applying the 
autocalibration technique (green dashed) that allows to access 
the systematic coefficients (gains and couplings for each chan- 
nel). 



each bolometer in the focal plane, the synthesized beam can be 
calculated very accurately. It can be compared to a direct mea- 
surement obtained by opening all the switches and scanning a 
source as well as to direct calculations using the positions of 
the horn and the optical simulation of the beam combiner. All 
these comparisons will allow for a number of consistency checks 
and measurements of systematic effects that will allow handling 
these effects in a very precise way. 

We have simulated this autocalibration technique in order to 
check its behavior (in the absence of noise) and we obtain ex- 
cellent results as shown in Fig. 10 We start from a set of visi- 



bilities observed with a 144 horns bolometric interferometer and 
apply randomly drawn systematic effects on each channel (com- 
plex gains and coupling using a Jones matrix formalism), the 
reconstructed visibilities are clearly corrupted (red dashed line) 
without using the autocalibration technique, but are exactly re- 
constructed (there is no noise in this simulation) when one in- 
cludes in the visibility reconstruction the systematics coefficients 
determined by the autocalibration technique (green dashed line, 
superimposed to the black one). This shows that the autocalibra- 
tion techniques actually works fine and allows recovery of the 
systematics in an efficient manner. An over- simplistic back of 
the envelope calculation shows that with NET - 300 yuK.sec 1/2 
detectors and a T ~ 100 K polarized source, one can reach an 
accuracy of order T/NET ~ 3 x 10" 6 on each of the Jones matrix 
coefficient if one spends one second on each baseline. For the 
whole array this would imply an autocalibration procedure that 
would last ~ 9 hours. It could be done once in a while and on 
a more regular basis one could perform the autocalibration only 
with baselines within subarrays of a smaller number of horns 
allowing to perform the autocalibration in about one hour. 

There is no equivalent to this autocalibration technique with 
an imager that will just rely on scanning sources (polarized and 
unpolarized) with the instrument. We can gather more informa- 
tions than that with our bolometric interferometer and therefore 
hope that the bolometric interferometer will allow to handle the 
systematic effects in a more accurate way than an imager, mak- 
ing the best of its interferometric nature. 
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7. Conclusions 

We have presented in this article the QUBIC instrument, a bolo- 
metric interferometer combining the signals from an array of 
wide (~ 14 degrees) entry horns using a cold telescope as an op- 
tical combiner to form interference fringes of /, Q and U Stokes 
parameters (modulated with a half-wave plate) on a bolometer 
array. We have shown that such an instrument can be achieved 
with technology that is mostly already available and that a small 
amount of R&D is required for some of the components, but 
that are common with future imagers. The image we observe on 
the focal plane is actually the so-called "synthetic image" mak- 
ing our bolometric interferometer a synthetic imager in contrast 
to usual interferometers whose observables are the visibilities 
(Fourier Transform of the synthetic image). With this synthetic 
imager the usual imaging techniques can be applied easily such 
as scanning the sky (in order to improve the sky coverage uni- 
formity and modulate /, Q and U), map-making and classical 
power- spectrum estimation. This simplifies the analysis with re- 
spect to the visibility-space one that is more difficult to achieve 
outside the flat-sky approximation. The statistical sensitivity of 
QUBIC has been shown to be comparable to that of an imager 
with the same number of horns covering the same sky fraction 
while we expect the systematic effects to be different from those 
of an imager. We also propose an autocalibration technique, spe- 
cific to bolometric interferometry, inspired by the redundancy 
based traditional interferometry techniques, that will allow to 
determine accurately the systematic effects (as models as Jones 
matrices) for each of the channels of the instrument. We antic- 
ipate that this will allow QUBIC to achieve a good control of 
systematics and therefore approach its statistical sensitivity. The 
full QUBIC instrument will comprise three frequencies and tar- 
gets to constrain at the 90% confidence level a tensor-to-scalar 
ratio of 0.01 with one year of data taking from the Concordia 
Station at Dome C, Antarctica. 
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